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Abstract 

We compute the decay constants of the pseudoscalar mesons D and 
D s using a linear combination of finite energy sum rules which min- 
imize the contribution of the unknown continuum spectral function. 
We employ the recent three loop calculation of the pseudoscalar two- 
point function expanded in powers of the running charm quark mass. 
The theoretical uncertainties arising from the QCD asymptotic ex- 
pansion are quite relevant in this case due to the relative small scale 
of the charm mass. We obtain the following results: fr> = 177 ± 21 
MeV and f Ds = 205 ± 22 MeV. These results, within the error 
bars, are in good agreement with estimates obtained using Borcl 
transform QCD sum rules, but somewhat smaller than results of 
recent lattice computations. 
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1 Introduction 



The decay constant of a pseudoscalar meson D q consisting of a heavy c-quark 
and a light quark q — u,d, s is defined through the matrix element of the 
corresponding pseudoscalar current as follows: 

<Q\(m c + m q )qi-f 5 c(0)\D q >= f Dq M% q . 

Recently, two different experimental groups have extracted the decay con- 
stant from a direct measurement of the absolute branching fraction for 
the Cabbing-suppressed leptonic decay mode D + — ► l + v ^E]- The results are 
quite different, although there is some overlap of the large experimental errors. 
Both results are compatible with the upper limit of 290 MeV, established by 
the MARK-III collaboration at 90% C.L. For the Df meson, there are 
several measurements of its decay constant fn+ published in the last decade 
H El El □ El M, with results in the range of 194 - 430 MeV. On the theo- 
retical side, there is a recent estimate based on Borel transformed sum rules 
fn = 195 ± 20 MeV [TU], and a preliminary result fjj = 225 + 11 — 13 i 21 
MeV ^1] obtained in three flavor lattice QCD. There appears to be some room 
for improvement of the accuracy of these results as well as for the systematic 
study of their uncertainties. On the other hand, the result /dJId = 1-18 for 
the ratio of the decay constants appears to be well established in lattice QCD 

[TT1II2 

In this letter, we estimate the decay constants fn and Jd s of the peudoscalar 
charmed mesons using an alternative method, based on finite energy sum rules, 
which compare moments of available experimental data with the corresponding 
QCD theoretical counterpart. In particular, we take linear combinations of 
positive moments in such a way that the contribution of the data in the region 
above the resonances turns out to be practically negligible. On the theoretical 
side we use a large momentum expansion of massive QCD at three loops. This 
expansion is known up to the seventh power of rn^/s, where m c is the mass of the 
charm quark and s is the square of the CM energy |13| . The expansion makes 
sense as long as s is far enough above the continuum threshold and above the 
resonances. On the phcnomcnological side of the sum rule we consider only the 
lowest lying pseudoscalar _D g -meson, once the unknown continuum contribution 
has been removed by a judicious use of quark-hadron duality in our method 

El 

The plan of this note is the following. In section 2 we briefly review the 
finite energy sum rule method employed. In section 3 we discuss the theoretical 
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and experimental inputs and we present our estimates for the decay constants. 
Finally, in section 4 we write up the conclusions. 



2 The method 

The two point function associated with the pseudoscalar current is: 

IL(s = q 2 ) = i ( dxe^ < n\T(j 5 (x)j 5 (0)) \n >, (1) 



where < | is the physical vacuum and the current j$ (x) is the divergence of 
the axial- vector current: 

j 5 (x) = (Mq + m q ) : q(x) i 75 Q(x) : (2) 

Mq is the mass of the heavy charm quark Q(x) and m q is the mass of the 
light quarks q(x), up, down or strange. The starting point of our sum rules is 
Cauchy's theorem applied to the two-point correlation function n(s), weighted 
with a polynomial P(s): 

/ P(s) II(s) ds = 0. (3) 
2tt% J r 

The integration contour T extends over a circle of radius \s\ = sq, and along 
both sides of the physical cut [s p /i ys .,So] where s p hy S . is the physical threshold 
and So is a duality parameter to be fixed with some stability criteria. The 
polynomial P(s) does not change the analytical properties of n(s). Therefore 
we obtain the following sum rule: 

- [ P(s) ImII(s) ds = — / P{s) n(s) ds. (4) 

n Js phys . 2m J\ S \=S0 

On the left hand side of equation Q we use the experimental information 
for Imn DATA (s) between the physical threshold s p h ys . and sq, whereas on the 
right hand side we use the asymptotic expansion n^ CD (s) of QCD, including 
perturbative and non-perturbative terms. The QCD expansion constitutes a 
good approximation of the two-point correlate on the circle for a large enough 
integration radius s . 

The experimental data are dominated by the first pseudoscalar cq resonance. 
In the narrow width approximation, the absorptive part of the two-point func- 
tion Imn DATA (s) can be split in two terms, the contribution of the resonance 
and the contribution of the hadronic continuum ImII cont starting at the physical 
continuum threshold s C ont. > Sphys. = M| , as follows: 

- Imn DATA ( S ) - M% fl 5{s - M 2 D ) + - Imff™' 6(s - W ) (5) 

TT q q q TT 

where Afo q and /u q are respectively the mass and the decay constants of the 
lowest lying pseudoscalar meson D q . 
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For the QCD correlate we write the decomposition, 

n QCD (s) = n pcrt -(s) + n nonport -(s), 



(6) 



We employ the two-point correlation function IT pcrt (s) with one massless and 
one heavy quark given to second order (three loops) in the strong coupling 
constant a s and expanded in a power series in the pole mass of the heavy quark 
including terms of order (M^/s) 7 . In the following compact expansion of 
the two-point function in terms of the pole mass M c can be found 



n port -(s) = n(°)(s) + 



a s (M c ) 



n«( s ) 



a s (M c ) 



n (2) (s), ( 7 ) 



where the different terms of the expansion in a s have the form: 



-1 k=0 



2 \ J 



(8) 



In the equations J7JSJ, M c is the pole mass of the charm quark. The coefficients 
are explicitly given in For instance, the one-loop term of the expansion 
in a s reads : 



n(°)( s ) = 



(M c + m q y S < 3-21og 



16tt 2 
3 + 2 log 

i /Mr r> 



M 2 
-4^ log 
s 



^ M 2 ) 



— s 

w 



m?\ 2 2 /Mr 3 



15 



1 f M. 



30 



s 

2 x 6 



3 V s 
2 (Ml_ 
s 



1 Mi 



6 



105 



The non-perturbative terms in the asymptotic expansion of equation © are 
due to the quark and gluon condensates. We will include terms up to dimension 
six [TSUinj: 



n „onpert. (s) = ^ + )2 ^ 



7T -S + M 2 



1 



1 



12 s- M, 



f2 



—G 2 



27 



-M c 
2 

M 2 



1 



Mt 



(s-Mf) 2 (s-M c 2 ) 3 
M 4 

u s (qq)' 



(qaGq) 



{s-M 2 c ) 2 {s-M*Y {s-M 2 c y 



For the quark condensate we include the a s correction |29|. it turns out to be 
small but non-negligible. 

In order to improve the convergence of the perturbative expansion, we re- 
place the pole mass by the running mass using the 0{a 2 ) result relating the two 
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El El El EDI- The perturbative piece of order (as)* of equation © can be 
rewritten in the form 



j=~l k=Q 



(9) 

and similarly for the non-perturbative piece. The coefficients Aj k depend on 
the mass logarithms ln(m^//i 2 ) up to the third power. As n(s)^ CD is not known 
to all orders in a s , the results of our analysis will depend to some extend on 
the choice of the renormalization point /i. In the sum rule considered here there 
are two obvious choices, fi = m c and // = t/sq. The former choice will sum 
up the mass logs of the form ln(m^//i 2 ) and the latter choice the ln(— s//i 2 ) 
terms. For definiteness, we take /i = m c . With this, the convergence of the 
perturbative terms is reasonably good. The results differ from taking \i = ^/sq 
by an amount consistent with the general three-loop asymptotic uncertainties, 
as we will analyze below. 

Looking back to equation |0J and taking all the theoretical parameters as well 
as the mass of the D g -meson and the physical continuum threshold as inputs of 
the calculation, we see that the decay constants can be computed from equation 
(JSJl only if we have good control over the hadron continuum contribution of the 
experimental side. 

To cope with this problem we make use of the freedom of choosing the 
polynomial in equation (0J. We take for P(s) a polynomial of the form: 

P n (s) = a + ais + a 2 s 2 + a 3 s 3 + . . . + a n s n , (10) 

where the coefficients are fixed by imposing a normalization condition at thresh- 
old 

Pn (s phys . = M 2 Dq ) = 1, (11) 

and requiring that the polynomial P n {s) should minimize the contribution of 
the continuum in the range [s con t., sq] in a least square sense, i.e., 

s k P„(s) ds = for k = 0,...n-l, (12) 

cont. 

The polynomials obtained in this way are closely related to the Legendre poly- 
nomials. In the appendix the explicit form of the set of polynomials used in this 
work is given. 

This way of introducing the polynomial weight in the sum rule minimizes 
the continuum contribution — Imn cont - 8{s — s cont .) on the phenomenological 
side of the sum rule. To the extend that Im n cont - can be approximated by 
an n-th degree polynomial these conditions lead actually to an exact cancella- 
tion of the continuum contribution to the left hand side of equation . The 
role of the D q resonance will be enhanced. We will see in our analysis that 
this choice of the polynomial has the additional effect of increasing the region 
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of duality characterized by the value of the duality parameter sq. In this way 
the asymptotic expansion of QCD can be used more safely on the circular in- 
tegration contour. Notice however that increasing the degree of the polynomial 
P n {s) will require the knowledge of further terms in the mass expansion and in 
the non-perturbative series. Therefore the polynomial degree cannot be chosen 
arbitrarily high. 

To check the consistency of the method, we have employed polynomials 
ranging from second degree to fifth degree, verifying that the results are com- 
patible within the range of the errors introduced by the incomplete knowledge 
of the QCD correlate and other inputs of the calculation. We also have checked 
explicitly that a smooth continuum contribution had no influence on the result. 

Our approach to suppress the continuum has been tested previously in the 
calculation of the heavy quark masses using analogue sum rules for the vec- 
tor current correlate where there exists more experimental information on the 
continuum. In the calculation of the charm quark mass, using the same poly- 
nomial method, the continuum, known from the BES II collaboration |21| . was 
shown to have practically no influence on the predicted mass |23| . Employing 
the same technique, a very accurate prediction of the bottom quark mass was 
also obtained using the experimental information of the upsilon system |24| . 

After these general considerations we proceed with the analytical calculation. 
The integrals that we have to evaluate on the right hand side of the sum rule, 
equation @, are 

J( P ,k) = ±J^(ln^) k dS) (13) 

for k = 0, 1, 2, 3 and p — —6, —5, .., n + 1. These integrals can be found e.g. in 
reference |23 ■ After integration, equation Q yields the sum rule 

M A Dq f D P{Ml q ) = m 2 M(m c {p) + m q ^)) 2 (14) 

p=0 i=0 j=-l fc=0 ^ ' 

+ non-perturbative terms 

where, for brevity, we have not written down the non-perturbative terms explic- 
itly. The contribution of the continuum is neglected as explained above although 
the continuum threshold is considered in the determination of the coefficients, 
a p of the polynomials (|12fl . 

Plugging the theoretical and experimental inputs (physical threshold, quarks 
and meson masses, condensates and strong coupling constant) into the sum rule, 
we obtain the decay constant fo q for various values of the degree n of the poly- 
nomial and various values of sq. Given the correct QCD asymptotic correlate 
and the correct hadron continuum, the calculation of the decay constant should, 
of course, not depend either on sq or on the degree n of the polynomial in the 
sum rule I0J. Accordingly, for a given n we choose the flattest region in the 
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curve /d„(so) to extract our prediction for the decay constant. To be specific 
we choose the point of minimal slope. On the other hand, for different polyno- 
mials, the value of fu , extracted in this way, could differ from each other as the 
cancellation of the continuum may be incomplete or the QCD expansion may 
not be accurate enough. We find, however, practically the same results for all 
our polynomials. This additional stability is truly remarkable as the coefficients 
of the polynomials of order n = 2,3,4 and 5 are completely different and the 
respective predictions are based on completely different superpositions of finite 
energy moment sum rules. This extended consistency leads us to attach great 
confidence in our numbers and associated errors. 

3 Results 

We calculate the decay constants for the D and D s heavy mesons. In the first 
case we take m q = everywhere. In the second case we retain m q — m s ^= in 
the factor (M c + m q ) 2 in front of the correlation function only. Further terms 
in the power series in m 2 s / s in (jSJ are completely negligible for the integration 
radii so we use in the calculation. 

The experimental and theoretical inputs are as follows. The physical thresh- 
old Sphys. is the squared mass of the lowest lying resonance in the cq channel. 
For q being the light quark u, we have: 

sphys. = M% = 3.493 GeV 2 (15) 

whereas the continuum threshold s cont is taken from the next intermediate state 
Dmr in an s-wave I = h, L e. 

Scont. = (M D + 2m,) 2 = 4.575 GeV 2 . 

For q being the strange quark we take: 

Sphys = A/ 2 s = 3.873 GeV 2 (16) 

The continuum threshold starts in this case at the value: 

Scont. = (M Ds + 2m x ) 2 = 5.009 GeV 2 . 

On the theoretical side of the sum rule we take the following inputs. The 
strong coupling constant at the scale of the electroweak Z boson mass |2*7] 

a s (M z ) =0.118 ±0.003 (17) 

that is run down to the renormalization scale using the four loop formulas of 
reference For the quark and gluon condensates (see for example an d 
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the mass of the strange quark [2S| we take: 

< qq > (2GeV) = (-267 ± 17 MeV) 3 , 

< — GG > = 0.024 ± 0.012 GeV 4 , 

< qaGq >= ml <qq>, with m 2 , = 0.8 ± 0.2 GeV, 
m s (2GeV) = 120 ±50 MeV, 

< ss > = (0.8 ±0.3) <qq > . (18) 

As discussed above, we fix the renormalization scale to be fi = m c {m c ). 
We use a reasonable variation of /i to analyze the corresponding uncertainty 
in our final result. Finally, for the charm quark, we take the value m c (m c ) = 
1.25 ± 0.10 GeV obtained by similar techniques |281 124| which is in a generally 
accepted range. 

In order to calculate the decay constant for the pseudoscalar meson D, we 
proceed in the way described above. We compute fjj as a function of so with the 
four different sum rules (@J corresponding to n = 2, 3, 4, 5. The results, plotted 
in Fig. 1 show remarkable stability properties. We define the optimal value of 
so as the center of the stability region (represented by a cross in Fig.l) where 
the first and/or second derivative of /d(so) vanishes. At these values of so we 
obtain the following consistent results: 

f D = 176 MeV for n = 2 (19) 
f D = m MeV for n = 3,4,5 

Notice from Fig. 1 that for the fifth degree polynomial (n=5) there is a 
stability region of about 20 GeV 2 around So , where the decay constant changes 
by less than three percent. From this change we estimate a conservative error 
inherent to the method of ±3 MeV. 

Other sources of errors arising in the calculation of fu are the quark con- 
densates which affect the result by ±5 MeV, and the charm mass which, in the 
range given above, produces a variation in the decay constant of —19, +9 MeV. 
This is one of the main source of uncertainty in the final result (see table 1). 
Considering the following results of the perturbative a s expansion: 

one-loop calculation: fffl — 142 MeV 
two-loop calculation: = 162 MeV (20) 

three-loop calculation: = 177 MeV 

we also take as a source of uncertainty the contribution of order a 2 , which 
amounts ten percent of the result. This yields an asymptotic error of ±15 MeV 
We point out that the convergence of the asymptotic series in the present 
calculation of the decay constant of the D meson is worse than the one we 
found for the B meson [55]. Finally we have considered the dependence on the 
renormalization scale in the range [i S [2,6] GeV. The error associated to this 
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f n (GeV) 




1U L\j JU H-U j 

S (GeV ) 

Figure 1: Decay constant fr> as a function of the integration radius sq for 
m c {m c ) = 1.25 GeV. The lines represent the sum rule (0} for several choices of 
the polynomial P n {s). 



change in fi is roughly related to the convergence of the asymptotic expansion 
and therefore it is not considered as an additional one. Other errors due to the 
QCD side of the sum rule, higher order terms inm^/s and the error on a s (mz), 
are negligible in comparison. 

From this analysis of errors, we finally quote the following result for the 
decay constant of D-meson 

f D = 177 ± 14(inp.) ± 15(asymp.) ± 3 (meth.) MeV. (21) 

The first error comes from the inputs of the computation, the second to the 
truncated QCD theory whereas the last one is due to the method itself. 

Proceeding in the same fashion, but keeping the mass of the strange quark 
in the overall factor and the order m s /m c in the one loop contribution, we find 
the decay constant for the D s meson, 



one-loop calculation: 


AO) 

JD B 


= 163 


MeV 


two-loop calculation: 


f (l) 
JD B 


= 188 


MeV 


three-loop calculation: 


f (2) 


= 205 


MeV 
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Inputs 


f D = 177 MeV 


m c (m c ) = 1.25 ±0.10 GcV 

< qq > (2 GcV) = (-267 ± 17 MeV) 3 
<^GG>= 0.024 ± 0.012 GcV 4 

< -qaGq > = ( 0.8 ± 0.2 GcV) 2 < qq > 
a s (m c ) = 0.335 ±0.010 


-19 ±9 MeV 
± 5MeV 
±lMeV 
±lMeV 
± 3 MeV 



Table 1: Uncertainties on fo from the uncertainties on the theoretical parame- 
ters 

and including the analysis of uncertainties we find: 

f Ds = 205 ± 14(inp.) ± 17(asymp.) ± 3(meth.)) MeV. 
(see Fig. 2) 

In the analysis of theoretical errors the only new ingredient is the uncertainty 
coming from the strange quark mass which turns out to be negligible. 

The ratio of the decay constants Jd 3 and fo (which would be 1 in the chiral 
limit) is of special interest. We find: 

^ = 1.16 ± 0.01 (inp.) ± 0.03 (meth.). (23) 

JD 

in complete agreement with lattice calculations. The uncertainties due to the 
theoretical inputs are correlated, so that the final error is very small. 

4 Conclusions 

In this note we have computed the decay constant of D g -mesons for q being 
either the strange or the u or d massless quarks. We have employed a judicious 
combination of moments in QCD finite energy sum rules in order to minimize 
the shortcomings of the available experimental data. On the theoretical side of 
the pseudoscalar two-point function, we have used in the perturbative piece an 
expansion up to three loops in the strong coupling constant and up to order 
(m^/s) 7 in the mass expansion and in the non-pcrturbative piece we considered 
condensates up to dimension six including the a s correction in the leading term. 
Instead of the commonly adopted pole mass of the bottom quark, we use the 
running mass to improve convergence of the perturbative series. 

In the sum rule, the contour integration of the asymptotic part is performed 
analytically. This particular fact differs from other computations based on sum 
rules where the asymptotic QCD is integrated along a cut of the two-point func- 
tion starting at the pole mass squared. The latter way to proceed is problematic 
when loop corrections are included and the complete analytical QCD expres- 
sion along the cut is not known. In this approach QCD has to be extrapolated 
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Figure 2: Decay constant fr> s as a function of the integration radius So for 
m c {m c ) = 1.25 GeV. The lines represent the sum rule for several choices of 
the polynomial P n {s). 



from low energy to high energy |29| . We also differ from many other sum rule 
calculation in that we do not require two unrelated sum rules to determine a 
duality point via an intercept of the curves /d(so). 

Our results are very sensitive to the value of the running mass, giving most 
of the theoretical uncertainty. They also turn out to be sensitive to variations 
of the renormalization scale /i. The uncertainties of other theoretical parame- 
ters like quark condensates and coupling constant are less important. Adding 
quadratically the different estimated errors we have the final results 

f D = 177 ± 21 MeV 

f D , = 205 ± 22 MeV (24) 

Comparing (|24(l with other results in the literature, our results agree within 
the error bars with the ones obtained using sum rule methods However, 
compared with lattice computations [341 135| they are a bit lower. 
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Appendix 



For convenience of the reader we list in this appendix the first few polynomials 
emerging from relations (|11I12() . From the second condition, namely 1|12|) . it is 
easy to realize that the set of polynomials P n (s) are n-degree orthogonal poly- 
nomials in the interval of the variable s € [s con t., so]- With the normalization 
condition Ijlll) (adopted to emphasize the contribution of the lowest lying reso- 
nance in the sum rule) they are related to the Legendre polynomials V n (x) in 
the interval of the variable x € [—1,1] as follows: 

Pn(s) = V f {3) \ (25) 



Where the variable x(s) is: 

x(s) 



2s - (Sp + Scont.) 
^0 -^cont. 



Obviously x(s) E [— 1,1] when s € [s con t.,so]- 

The explicit form of these polynomials is well known and can be found, for 
instance, in [3^]. Nevertheless, for sake of completeness, we quote here the ones 
we have used in the calculation. 

?2 (*(*)) = ^ 2 -l), 
V 3 (x(s)) = ~(5z 3 -3x), 

Va{x{s)) = i(35a; 4 -302: 2 + 3), 
8 

V 5 (x(s)) = l(63x 5 - 70x 3 + 15a;). 
8 

Finally, in Fig. 3 and in order to appreciate the suppression of the exper- 
imental physical continuum data in the sum rule, we have plotted the form of 
the polynomials P n (s) for n=2,3,4,5 at the stability values of s$ used in the 
calculation of fn- 
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